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Abstract. An abstract should be given 



The recent influx of high-quality infrared spectro- 
scopic data has prompted an extensive reassessment of 
various laboratory models in comparison with the ob- 
served UIBs (Unidentified Infrared Bands). As a result, 
significant modifications were brought to the original 
paradigms. The focus here is on the evolution of the coal 
model, characterized by 1) a shift towards less aromatic 
materials of the same family (kerogens), 2) the introduc- 
tion of a new excitation mechanism (chemiluminescence), 
based on the capture of hydrogen atoms by carbonaceous 
dust. Both developments are intended to accomodate ob- 
servations from a larger range of dust environments and 
evolutionary stages. This leads to a more quantitative de- 
scription of dust composition and structure, and a better 
understanding of its history. In short, according to the 
present model, the evolution of dust from inception in the 
circumstellar shells of AGB stars, through strong inter- 
stellar radiation fields, to consumption in protostars, is 
approximately mimicked by progressively more aromatic 
materials, starting from the young and mostly aliphatic 
kerogens, through more and more mature coals, to the fi- 
nal stage of polycrystalline graphite. A similar family of 
materials is obtained in the laboratory by annealing in 
vacuum up to or beyond 3000 K. The composition, struc- 
ture and IR spectrum of these materials are extensively 
documented in the geology literature. 

The fundamental characters of this model, viz. chem- 
ical and structural disorder, and diversity of chemical 
bondings, naturally point to ways of further tailoring, in 
order to fit particular observations more closely. 

Key words. Interstellar dust-UIBs-PAHs-coal/kerogen- 
chemiluminescence-IR spectra 



Send offprint requests to: R. Papoular 



1. Introduction 

This paper is devoted to the understanding of the set of 
celestial IR (InfraRed) emission bands in the range 3- 
15 /jm which have been loosely, but wittingly, dubbed 
Unidentified Infrared Bands (UIBs) well after their dis- 
covery in the '70s. Since then they have been thor- 
oughly documented, using higher sensitivity and resolu- 
tion. Particularly large amounts of high quality data were 
collected by the IR satellites ISO [1] and IRTS [2], which 
were launched and operated in the '90s. Knowledgeable 
reviews of these and other observations may be found, for 
instance, in [3,4,5]. 

While electronic transitions in atomic ions, such as 
Argon and Sulphur, are also detected in the relevant range, 
sometimes on top of a UIB, they are easily distinguished 
by their much narrower width. The only present consensus 
regarding the UIBs is that they are due to fluctuations of 
the electric dipole moment associated with the vibrations 
of the atoms of carbon-rich carriers. Different research 
groups do not agree as to whether the carriers are free- 
flying molecules or solid-state grains. Very closely related 
to this dichotomy is the issue of the excitation mechanism 
which gives rise to the IR emission. Roughly speaking, 
the free-flying molecule model has been associated from 
inception to stochastic heating, i.e. random absorption of 
a single UV photon, followed by immediate thermaliza- 
tion of the deposited energy, and subsequent cooling due 
to IR radiation; this of course requires very small carriers 
(<50 atoms) and strong UV irradiation. Moreover, these 
carriers are considered to be PAHs (Polycyclic Aromatic 
Hydrocarbons), which are planar clusters of hexagonal C 
rings with peripheral attached H atoms. No assignment 
to any known particular PAH has yet been proposed, but 
it is suggested that the problem will ultimately be solved 
by a size-distribution of neutral/ionized, more or less hy- 
drogenated PAHs. In this model, separate particles, the 
Very Small Grains (VSG), have to be invoked in order to 
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account for the observed underlying continuum, but no 
composition or structure has yet been specified for this 
family as well. Donn et al. [6] have discussed this model 
long ago, but Sellgrcn, Tokunaga, Boulanger and others 
[5,3,4] have recently spelled out some of the "puzzles" [5] 
presented by the new observations. 

In this paper, I explore the extent to which these puz- 
zles can be solved by a solid-state model [7], and more 
specifically the coal/kerogen model. This differs from the 
PAH model on nearly every count. The composition, here, 
is not restricted to C and H, but also includes small 
amounts of O, N and S, as suggested by cosmic abun- 
dances. The structure is not planar or regular, but 3-D 
and highly disordered; it is partly aromatic, partly aliphatic 
and olefinic depending on dust environment and history. 
To accomodate this diversity and disorder, the carrier size 
should be larger than a few hundred atoms. Most im- 
portantly, well defined natural terrestrial analog materials 
are proposed: these are all members of the kerogen/coal 
"saga", which have been thoroughly documented during 
the past half century, and whose composition and struc- 
ture have been abundantly described in the literature. The 
diversity of the celestial IR spectra, in relative feature in- 
tensities and in position and width, is matched by that 
of the terrestrial spectra, which itself is not erratic but 
known to be due to progressive release of S, N, O and 
H atoms upon ageing and/or heating of the material. No 
ionization or size distribution need be invoked. The struc- 
ture accounts for the feature bandwiths and ensures the 
presence of a continuum where required. Finally, the exci- 
tation mechanisms invoked here are 1) steady state "ther- 
mal" heating for strongly illuminated dust, 2) chemilumi- 
nescence for weakly illuminated dust. 

The coal model was first presented in 1989 [8]. At that 
time, modelling the available celestial spectra required the 
use of relatively highly aromatic coal analogs [9] . The data 
now available warrant the complementary use of the lesser 
aromatic analogs, kerogens [10]. These data also prompted 
the consideration and elaboration of a new excitation pro- 
cess [11]. I therefore believe that the coal/kerogen model 
is now in a position to cure the main riddles of the PAH 
model. These developments motivated the writing of the 
present report. 

Section 2 is a short reminder of the main properties and 
behaviours of kerogens and coals, relevant to the present 
discussion. In Sec. 3, I delineate some of the main riddles 
signalled by authors who applied the PAH model to the 
interpretation of astronomical observations. In each case, 
I then discuss how the solid-state model can remedy the 
problem or provide an alternative interpretation. Section 
4 discusses in more detail chemiluminescence as an alter- 
native to stochastic heating. Operational conclusions are 
drawn in the last section. 

2. Some properties of coals and kerogens 

The following is a cursory summary of the properties of 
natural amorphous carbons of interest to the present issue. 
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Fig. 1. Van Krevelen diagram for coals and kerogens. The 
strip defined by the two continuous lines is the locus of 
coals of various origins and depths; I, II and III refer 
to kerogens of lake, marine and continental origins, re- 
spectively; the arrow shows the initial general direction of 
spontaneous changes in H and O atomic concentrations; 
see text and [7] for details. 



Details and abundant references to specialized literature 
can be found in [7,10]. 



2.1 Composition and Structure 

The term coal refers here to the residual part of the 
raw terrestrial material after minerals and free molecules 
preexisting in its pores have been disposed of by suitable, 
standard, physical and chemical procedures. This residue 
is almost exclusively composed of C, H and O, with traces 
of N and S. While coal is extracted from deep in the earth, 
kcrogen is obtained from the superficial sedimentary rocks 
after treatment by aqueous alkaline and other organic 
solvents. The solid residue can be considered as coal in 
"dispersed" form (very small particles). It is in this form 
that most of the organic matter in meteoritic carbonaceous 
chondrites is found, a useful hint for our present purposes. 

The properties of coal do not depend heavily on the 
geographical origin of the material. As a result, each sam- 
ple is concisely characterized by one point in the van 
Krevelen diagram (Fig. 1) with coordinates O/C and H/C 
(atomic ratios). The representative points of coals from 
different origins form an inverted- L-shaped strip. For coals 
extracted from various mining depths, but same geograph- 
ical location, the strip nearly reduces to a line of the same 
shape; as the mining depth increases, the representative 
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point follows the line down to the origin. The composition 
and properties of kcrogens are more dispersed at high het- 
eroatom concentrations, but ultimately merge with those 
of coals at low heteroatom concentrations. This graphical 
orderliness stems from the chemical properties of carbon. 

Like in PAHs, the atoms in coals and kerogens are 
held together by a and 7r covalent bonds. But, while the 
former include only aromatic C sites, the latter comprise 
all three types of sites: sp 1 (acetylenic) , sp 2 (aromatic, 
olefinic) and sp 3 (aliphatic). This is a result of the pres- 
ence of "impurities" (O, N, ...) which impedes the tcn- 
tency towards the more thermodynamically stable aro- 
matic structures and leads, instead, to disordered materi- 
als. As a consequence, PAHs know of only one type of C-H 
bond, while coals/kerogens admit of several others: C-H 
(non-aromatic), CH 2 (olefinic, aliphatic) and CH 3 . Other 
important consequences are: a diversity of C-C bonds (sin- 
gle, double), the presence of COH and C=0 groups (which 
all have strong incidence on IR spectra), and -O- bridges 
which provide the structure with some flexibility and allow 
it to grow in 3-D. 

Together with disorder comes the potentiality for evo- 
lution. Figure 2 illustrates the displacements brought 
about, in van Krevelen's diagram, by spontaneous or in- 
duced losses of CO2, H 2 and CH 4 . In earth, these losses 
are stimulated by heat and pressure, which explains the 
trajectory of the representative point in Fig. 1 towards 
the origin as mining depth increases. The structural rear- 
rangements which necessarily accompany these composi- 
tional changes are schematically outlined, in Fig. 2, along 
the left vertical axis, showing the aromatization trend 
along the evolutionary (maturation) track. 

In the laboratory, mild, prolonged, heating (anneal- 
ing) of a low aromaticity sample prompts its represen- 
tative point to quickly proceed along the same natural 
track it would have followed in earth, given time, had it 
not been extracted. This simple and coherent behaviour 
of coal/kerogen again emphasizes the tight link between 
composition and structure and provides a guiding line to 
find one's way in the diversity of these materials. 

Thousands of samples have been studied thoroughly 
all over the world, using all available analytical tech- 
niques: thermal, thermogravimetric, IR, Raman and NMR 
spectroscopies, etc. This made it possible to find correla- 
tions between spectral bands and assign nearly every band 
to a functional group (a handful of associated atoms). 
These efforts culminated in the building up of represen- 
tations of the structures of coals and kerogens as a func- 
tion of their evolutionary stage, i.e. their location in the 
van Krevelen diagram (Fig. 1). The work of Behar and 
Vandenbroucke [12] is a good example of this develop- 
ment. They found that it is not possible to adequately de- 
scribe the observed spontaneous and continuous changes 
in properties in terms of a limited number of small, spe- 
cific, molecules (a conclusion to be extended below to the 
interpretation of UIBs). In order to fit the measurements 
on different samples, one has rather to build a random ar- 
ray made of a dominant carbon skeleton with functional 



groups of heteroatoms (H, O, N, S) attached randomly to 
it, then statistically tailor the number of different bond- 
ings (C-C, C=C, etc), functional groups (C-H, C=0, etc) 
and aromatic or polyaromatic rings, etc. The variety of 
environments of the functional groups ensures that their 
characteristic vibrations will blend into bands of the right 
width. The relative intensities of the latter will change ac- 
cording to the concentrations of the corresponding func- 
tional groups. This scheme accounts for the continuous set 
of IR spectra observed along the evolutionary track. 

Behar and Vandenbroucke give numbers for the struc- 
tural parameters enumerated above and for 8 represen- 
tative evolutionary stages. They also give sketches of the 
corresponding structures, examples of which are shown in 
Fig. 3. A remarkable feature in these sketches is that, even 
though the size and number of aromatic clusters are quite 
limited, all the main UIBs are present in the corresponding 
IR spectra, albeit with different intensities. The drawing 
of the clusters also highlights the 3-D nature of the struc- 
tures. Note the large number of short aliphatic chains in 
the lesser aromatic samples; the aromatic clusters become 
large and dominant in the late stages of maturation. 



2.2 The IR spectrum 

As a typical example, Fig. 4 (adapted from [13]) dis- 
plays the IR absorbance spectrum of a young kerogen 
of type II (H/C=1.32, O/C=0.104). The peak wavenum- 
ber (cm -1 ) and wavelength (/im), integrated intensity 
(cm/mg) and assignment of the bands are, respectively, 
(see [14]): 

1. -3400(^2.95); 46.2; OH stretch; due partly to 
chemically bonded OH groups and partly to adsorbed or 
trapped, not chemically bonded, H 2 molecules; the lat- 
ter are responsible for the long redward tail, through H- 
bondings with other parts of the skeleton; peak position 
and band profile change considerably with H 2 content 
(depending on evolutionary stage or heat treatment). This 
feature is clearly distinct from the water ice band, which 
peaks near 3.1 /jm and has a much steeper red wing (as- 
sociated to that observed in the sky towards young stellar 
objects). Adsorbed H 2 is probably absent from IS dust. 

2. 3060 (3.27); C-H aromatic or olefinic stretch, barely 
visible here, and only measurable for very small H/C and 
O/C ratios. 

3. 2920 (3.42); 48.2; blend of anti-symmetric and sym- 
metric CH 3 stretch at 2962 and 2872, asymmetric and 
symmetric CH2 stretch at 2926 and 2853; CH stretch 
at 2890 cm" 1 , respectively. Here, the non-aromatic C-H 
bonds reside mainly on non- aromatic structures, not at 
the periphery of PAHs like in the work of Wagner et al. 
[15]. 

4. 1710 (5.85); 14; C=0 (ketone) stretch. 

5. 1630 (6.15); 11.4; disputed (but perhaps concurrent) 
assignments to H 2 deformation, quinonic C=0 with H 
bond and C=C olefinic and aromatic stretch. 
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Fig. 2. Van Krevelen diagram (after Durand [14], pl22). Each arrow represents the direction of loss of a particular 
gas molecule in the course of evolution. A rough sketch of the main structural elements is given on the left side for a 
few evolutionary stages. 



6. 1455 (6.87); 6.1; asymmetric CH 2 and CH 3 defor- 
mation. 

7. 1375 (7.28); 1.2; symmetric CH 3 deformation. 

8. 1800 to 900 ( 5.5 to 11); 67.3; massif (underlying 
broad band) peaking near 8 pra, due to C...C and C-0 
stretch, C-H in-plane bend and OH deformation; in highly 
aromatic samples, there may be a contribution by elec- 
trons, in the form of a plasma resonance [16]. 

9. 930 to 700( 11 to 14); aromatic out-of-plane bend- 
ing, depending on the number of adjacent protons; weakly 
contrasted 3 or 4 peaks, barely visible here; intensity in- 
creases with mining depth but never exceeds an intensity 
of 5. 

10. An underlying continuum, roughly decreasing with 
frequency, is mostly visible in the near-IR, but also present 
farther to the red in all coals/kerogens; part of it is due to 
light scattering (in transmission measurements) and part 
of it to real absorption; it increases strongly with aro- 
maticity. 



The intensity, K, given here (third number) is related 
to the absorbance, a, by 

a(cm- 1 )=2.3 10 3 p K/Av, (1) 

where p is the material density in g cm~ 3 (^1) and Au, 
the bandwidth in cm -1 . K (cm/mg) is 10~ 3 times the 
intensity usually defined in the astrophysical literature; it 
is also related to the integrated cross-section per C atom, 
S, by 

^(cm^m-iCatonr 1 ) = 4.7 10- 20 (1 + + ±§g) K (2) 

« 4.7 to 7 10~ 20 K. 

Not unexpectedly, a strong relationship holds between 
the bands at 2920 and 1455 cm -1 for all kcrogen types and 
evolutionary stages: the ratio of their integrated intensities 
is ~8. Also, K(3.4 pm) increases linearly with H/C, from 
at <0.3 at least up to 80 cm/mg at H/C=1.3. 

It should be clear from the above that a single spec- 
trum like Fig. 4 cannot, by far, convey the full diversity 
of coal and kerogen spectra. Series of spectra in restricted 
wavelength ranges and corresponding to successive degrees 
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Fig. 3. Chemical representation of a type II kerogen at an early (b), and a later (c) stage of aromatization (adapted 
from Behar and Vandcnbroucke [12]). Following common practice, C-H bonds are not represented. Aliphatic carbon 
chains (alkanes) are shown as broken, undulating, lines. The aromatic clusters of benzenic rings are shaded. Various 
functional groups and oxygen bridges are labeled. Note the increase in number and size of aromatic clusters, and the 
decrease in length of the aliphatic chains, from (b) to (c). 
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Fig. 4. Typical IR absorbance spectrum of a type II kerogen, after Robin and Rouxhet [13]). The dashed areas are 
used to determine the band intensities, K (see eq. 1). A number of partial assignments of bands are enclosed in squares; 
see details in Sec. 3. 



of maturation may help for this purpose. An example is 
given in Fig. 6 for the 3 jim range; others can be found in 
[10, 14]. 

3. Puzzles and clues 

Here I discuss some of the problems encountered in ap- 
plying the PAH model to the most recent astronomical 
observations. In enunciating these, I shall try to remain 
faithful to the wording of the referenced papers; this task 
is made easier by the authors' candour. The recent review 
by Sellgren [5] was particularly useful in this respect. 



3.1 Carrier composition 

By definition, PAHs are pure hydrocarbons. However, 
the carriers of the UIBs are most likely formed in the 
winds of C-rich post AGB stars. In general, oxygen is not 



much less abundant than carbon in these winds. Nitrogen 
and sulphur are also present in notable quantities. It is 
therefore difficult to understand why these three elements 
should be excluded from the dust formed in such envi- 
ronments. All the more so since the less abundant silicon 
shows up in the SiC feature at 11.3 /im, and, most im- 
portantly since the ISM (InterStellar Medium) harbors 
many families of molecules containing these so-called het- 
eroatoms, especially O, mostly associated with C (see 
[17,18]). 

On the other hand, minority atoms, mostly O, arc 
present in coals and kcrogens in known concentrations, 
and these inclusions have important consequences: 1) they 
considerably enhance the dipole moment of some vibra- 
tional transitions involving C atoms, such as the 6.2 ^m 
stretch of the C=C bond; 2) they give rise to features ob- 
served in the sky in the mid-IR, such as the band at 1710 
cm -1 (5.85 /xm) and the massif at 7-9 /tm [10] , as well as in 
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the far IR, such as the 30 blend associated with COH 
wagging (see [19]); 3) minority atoms efficiently interfere 
with the tendency of C atoms to settle into an ordered 
aromatic structure; hence the amorphous 3-D structure of 
Fig. 3, which determines all the favourable properties of 
this model, as shown below. 



3.2 Feature bandwidths 

The UIBs are strikingly wider than the bands of any 
observed and identified molecule, the narrowest width be- 
ing about 15 cm -1 . Even at the highest resolution, no 
rotation-vibration decomposition is discernible, although 
such decomposition was clearly demonstrated, with the 
same spectrometers, for many identified molecules, such as 
H2 and H2O. Neither is it possible to assign these widths 
to lifetime broadening, since the typical width of 30 cm -1 
would correspond to a lifetime of about 0.1 ps, way shorter 
than the radiative lifetime (~0.1 s) and non-radiative life- 
times, ranging between 1 and 1000 ps. Thermal Doppler 
effects are also excluded for the masses and temperatures 
of particles of interest. Thus, although UIBs can be fit by 
one or more Lorentzians or Gaussians, these are too wide 
to carry physical meaning: the UIBs are definitely broader 
than laboratory PAH features in absorption (see [20] , [5] ) . 

In emission, highly excited molecules may have much 
broader features than in absorption because of anhar- 
monicity and spectral congestion. This is illustrated by 
Williams and Leone (21], who excited gas-phase naph- 
talene molecules (plane clusters of 2 benzenic rings, 18 
atoms) with short laser pulses of photons of energy 5.3 and 
6.6 eV. Monitoring the outgoing radiation at 3.3 (im, (C- 
H stretch), they found its width to decrease from ^100 to 
~30 cm -1 , while the peak shifted to the blue by 45 cm -1 , 
as the energy content of the molecule was reduced by on- 
going radiative emission. No such behaviour is observed 
in the sky. 

Being associated by birth with single-photon excita- 
tion (stochastic heating), the PAH model is therefore faced 
with a dilemma: either the molecule is so small, and the 
photon so energetic, that the "temperature" can reach 
about 1000 K and the emission is strong enough as to be 
detected, but then with unacceptably large widths and red 
shifts; or the molecule is large and the photon energy low 
so that no undesirable broadening and/or red shift occurs, 
but then with unacceptably weak intensities and narrow 
widths. Recent attempts to solve this dilemma relied on 
the strong but unsupported assumption that the band po- 
sition and widths of the features are the same functions 
of temperature for all interstellar PAH sizes, or on very 
specific mixtures of selected PAH cations supposed to be 
independent of history and environment (see discussions 
by Verstraete et al. [22, 23] and Boulanger et al. [24]). 

By contrast, the disordered solid-state model assigns 
the features to the synchronous vibrations of specific func- 
tional groups, and most of their width to the dispersion of 
similar functional groups among various environments, a 



situation corresponding to small-scale inhomogencity [20] . 
The widths of UIBs is indeed typical of solid, disordered 
materials [25] and, in particular, coals and kerogens [10]. 
The disorder, here, is mainly induced by heteroatoms and 
is, therefore, chemical and structural in nature (Sec. 2). 
Within limits, and especially below 500 K, feature widths 
are not sensitive to temperature or grain size, provided 
it is large enough (a few hundred atoms) to accomodate 
the required diversity of environments for a given average 
chemical composition. Under such circumstances, small 
degrees of ionization and dehydrogenation arc not likely 
to alter the spectrum notably, in stark contrast with PAH 
behaviour. 

While this explains the remarkable stability of the UIB 
widths [3,22], it also leaves room for the observed diver- 
sity [20], which may result from different compositions at 
inception, as well as different thermal and radiative pro- 
cessing in the interstellar medium. Because of this disor- 
dered, solid-state structure, the coal/kerogen model also 
adequately mimicks the celestial massifs at 7-9 /um (see 
fig. 4) and 11-13 fim [9], which are major puzzles for the 
PAH model [5]. 



3.3 Carrier aromaticity 

The use of the term PAH implies full aromaticity. This 
radical assumption was initially suggested by the strength 
and ubiquity of the 3.3, 6.2 and 11.3 /im features. However, 
it was not supported by later observations illustrated by 
fig. 5 for the stretch bands and in [7] for other ones: dis- 
tinct emission features between 3.3 and 3.6 /im are seen to- 
wards AGB stars such as IRAS 05341+0852, 04296+3429 
and CRL 2688, in Planetary Nebulae (e.g. Hb 5 [26]) and 
towards at least one Nova: Nova Cen 1986 [27]. In some 
cases, they are even stronger than the 3.3 /xm feature, as 
in absorption towards the Galactic Center, where the aro- 
matic feature is dwarfed by a massif with peaks near 3.4 
and 3.5 ^m. 

The initial assignment of these "red" features to hot 
bands (v>l) of the C-H str transition has proven uncon- 
vincing because of their strength and the non-detection of 
the correlative first harmonic to be expected near 1.67 /zm 
[15, 26]. It seems therefore more reasonable to assign them 
to the well documented aliphatic and olefinic, symmet- 
ric and asymmetric, C-H stretches at 3.4, 3.46, 3.515 and 
3.56 /mi, in CH, CH 2 , CH 3 and aldehydes (-CHO) groups 
[33,34]. Figure 6 illustrates this statement by setting, side 
by side with the spectra of fig. 5, the spectra of coals at 
increasingly higher stages of aromatization (see Sec. 2): 
not only is the fit of corresponding spectra remarkable, 
but the regular decrease of the red blend together with 
the rise of the 3.3 feature confirms their assignment to 
non-aromatic vibrations. 

Neither is it obvious that the 6.2 and 11.3 fim features 
are characteristic of aromatics. Indeed, while the 11-13 /im 
massif with 3 or 4 peaks, observed towards post- AGB stars 
may reasonably be ascribed to the out-of-plane vibrations 
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Fig. 5. The C-H stretch band of a) GC/IRS6E, in absorp- 
tion (Pendleton et al. [28]); b), c), d) 3 post-AGB stars 
in emission: IRAS 08341+0852 (Joblin et al. [29]), IRAS 
04296+3429 (Geballe et al. [30]) and CRL2688 (Geballe 
et al. [30]); e) reflexion nebula NGC 2023 in emission 
(Sellgren [31]). 



of aromatic C-H bonds [9], it has proven very difficult 
to devise a PAH cluster that could display a highly con- 
trasted, strong and lonely 11.3 feature as observed in most 
IS (InterStellar) spectra. According to [35], this requires 
"rectangular" PAHs with 100 to 150 C atoms, a very spe- 
cial requirement indeed! (see next subsection). Moreover, 
the expected correlation between the 3.3 and 11.3 fea- 
tures is not confirmed by observations [47, 52, 63]. Even 
the exact spatial coincidence of the intensity peak of the 
11.3 feature with that of the 3.3 feature appears to be in 
doubt (see [63, 64]). 

On the other hand, a (more familiar) alternative or 
concurrent assignment for the 11.3 (im band could be 
the wagging of the two H atoms of a vinylidine group 
(R2-C=C-H2), which also occurs at the desired frequency 




Fig. 6. The C-H stretch absorption band of coal as a func- 
tion of evolutionary stage: a) Vouters mine (O/C=0.06, 
H/C=0.75), b) Mcricourt mine (O/C=0.028, H/C=0.59), 
c) Escarpelle mine (O/C=0.018, H/C=0.46), d) and e) 
Escarpelle sample annealed at 525 and 600 K, respectively. 
Adapted from Guillois [32]). 



[33] . Clearly, the constraint of full aromaticity is not war- 
ranted: an understanding of the "3.4" feature, as well as 
of the diversity of UIB spectra, is better served by think- 
ing in terms of evolution of IS dust composition and struc- 
ture, both intimately connected and induced by prolonged 
heating (Sec. 2) and irradiation [36]. While this is natu- 
ral, straightforword and quantitative in the coal/kerogcn 
model, it seems to be hardly possible with compact PAHs, 
whatever their size and degree of ionization or hydrogena- 
tion. 



3.4 Spectral continuum 
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Known PAHs, even when ionized, cannot provide 
enough visible continuum absorption to account for the 
energy in the emission features observed around UV-poor 
reflection nebulae [5]. Neither can they mimick the unpo- 
larized near-IR continuum detected towards the same ob- 
jects [5,3]. This, of course, does not come as a surprise 
since free-flying molecules are not usually known for their 
continuum, although they do exhibit lines or bands in the 
vis/UV. Herein lies yet another, perhaps the main, rea- 
sons for invoking PAHs much larger than initially envi- 
sioned [5,35]. Unfortunately, none has yet been produced 
in the laboratory, if this is possible at all. Indeed it is well 
known that, in the absence of externally applied pressure, 
a planar sheet of compact benzene rings (graphene) of in- 
creasing size tends to bend, warp and finally folds upon 
itself to form a fullercne or a carbon nanotube[37]. 

Now, these particular carbon structures indeed have 
features in the UV/vis but they also have features else- 
where, such as the four conspicuous IR bands of fullcrene 
at 526, 576, 1183 and 1428 cm" 1 [38]. None of these has 
been detected in the sky despite intensive seach. 

Moreover, even these structures do not display mid- 
and far-IR continua, which points to still another puzzle 
of the PAH model: there is also a mid-IR continuum un- 
derneath the UIBs; for want of a better assignment, this 
has been attributed to a separate dust component, the 
very small grains (VSGs), but its composition and struc- 
ture could not be specified further [5, 39]. Note that they 
were thought to be small solely on the grounds of the self- 
imposed constraint that their emission should be excited 
upon absorption of a single UV/vis photon. Recent sur- 
veys of the spatial variations of this continuum have added 
to the confusion and seem to require again very special 
behaviour of the assumed separate component in order to 
match the systematic changes in the ratio of features to 
continuum [24, 39, 40, 41, 42]. 

Of course, theorists never lack an explanation, but all 
tentative applications of the PAH model seem to converge 
towards particles much larger than invoked in the first 
place [5, 24, 35, 43, 44]. At that point, it must be recalled 
that bulk carbonaceous materials indeed display a contin- 
uum, but the profile and intensity of the latter depend 
heavily on the composition and structure of the material 
[7]. Thus, the continuum of glassy carbon is very strong 
and varies approximately as A -1 at least as far as 100 /zm, 
but coals and a-C:H (amorphous hydrogenated carbon) 
are semi-conductors whose absorption falls down sharply 
towards the red at a cut-off wavelength (~1 /um) which 
increases with aromaticity. The lower aromaticity kerogen 
has a cut-off at shorter wavelengths and, hence, a weaker 
mid-IR continuum, better suited to the high contrast of 
the UIBs; it has also a near-IR continuum rising towards 
short wavelengths; it has the further advantage of carrying 
both continuum and bands [10], thus eschewing the need 
for two different carriers. 



3.5 Near- and far- UV bands 

Simple aromatic molecules are known for their near- 
and far-UV transitions between ground and excited elec- 
tronic states[45]. None of these was detected in the sky, 
either in absorption or in emission (fluorescence) [6], al- 
though several attempts were made to identify the DIBs 
(Diffuse Interstellar Bands) with such bands. In this field, 
too, hopes now hinge upon large 3-D systems [46]. 

Indeed, it was shown that high-temperature treatment 
(HTT) of coals favored aromatization and the growth 
of a single near-UV feature which irreversibly becomes 
narrower and bluer [7]. In some favourable cases, the 
end point of this evolution is poly crystalline (disordered) 
graphite, with a single, conspicuous, feature near 217 nm, 
on top of a continuum which increases sharply to the blue 
and decreases slowly to the red [7]. Then, of course, the 
material can hardly be designated as a PAH. 



3.6 Infrared bands 

Regarding wavelength and bandwidth, the UIBs are 
remarkably stable against large variations in the environ- 
ment and irradiation. This stability is incompatible with 
an assignment of the bands to a soup of free-flying inde- 
pendent PAHs, neutral or ionized and more or less hy- 
drogenated, as suggested by a number of authors in the 
past, since later laboratory studies showed that individual 
PAHs arc highly sensitive to irradiation wavelength and 
intensity (see discussions in [5, 24, 48, 49, 50]. In trying to 
use these parameters to fit observations, one is often led 
to conflicting or paradoxical conclusions (see discussions 
in [22, 23,35, 49, 50, 51]). 

One is therefore led again to prefer a solid-state (bulk) 
model: in that case, ionization and dehydrogenation can 
only affect the surface of the grain, a minor fraction of the 
whole. Also, the dust is then much more robust against 
destruction by shocks. Moreover, the structure and com- 
position of a bulk, amorphous, material as coal/kerogen 
can evolve slowly, spontaneously or upon heat treatment 
as shown in Sec. 2. This is accompanied by subtle changes 
in band wavelengths and widths (small fractions of a mi- 
crometer), also observed in the sky (see [10, 20]). Very 
large changes in relative intensities are also observed, espe- 
cially in the ratios [3.3]/[3.4], [3.3]/[11.3] and [11.3]/[12.7] 
(see[47, 52]). On terrestrial analogs, most of these can be 
explained by coherent and systematic changes of compo- 
sition and structure brought about by heat or slow spon- 
taneous evolution. The most spectacular and best docu- 
mented case is that of C-H stretch (see fig. 6): Dischler 
analyzed in detail the evolution of each component of this 
group of vibrations [34] . 

Thus the spectral diversity of coals /kerogens of differ- 
ent degrees of maturation ( ranks ) and origins matches the 
diversity of astronomical spectra. The main astronomi- 
cal features associated with carbonacious IS dust, includ- 
ing the continuum, are provided, to various extents, by 
coals /kerogens of various ranks. Different coals /kerogens 
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or mixtures thereof can account, roughly at least, for as- 
tronomical spectra from different IS environments. These 
models introduce no significant spurious feature anywhere 
in the explored spectrum. Examples of quantitative fits be- 
tween UIBs and models can be found in [7] for the UV/vis 
range, here in Fig. 5 and 6 for the near-IR, in [10, 11] for 
the 5-10 /im range, in [9] for the 11-13 ^m range and [19] 
for the far-IR range; see also [32]. 

That is not to mean that a perfect match between 
celestial and laboratory spectra is easily obtained. For one 
thing, the high quality of ISO spectra has revealed many 
new weak features, or details of the profiles of previously 
known features, which have yet to be elucidated (see [35, 
52]). More importantly, no satisfactory interpretation is 
yet available for the 12.7 feature, nor for the 11.3 feature 
being usually so much stronger than the underlying blend 
(11-13 /urn) which is commonly assigned to C-H mono, 
duo, trio and quarto out-of-plane bends (see Sec. 2). The 
11.3 and 12.7 bands are not present with a high contrast 
in the solid- state model spectra. 

But the coal/kerogen model provides at least a frame- 
work with definite guidelines and valuable assets: 

-representative natural samples are available with 
known composition, structure, physical properties and 
band assignments; 

-well-defined heat treatments can alter these properties 
controllably and coherently; 

-the range of possible model spectra can be extended 
by using samples of different depths and geographical ori- 
gins, thus scanning the whole van Krevelen diagram (fig. 
1), i.e. a large range of natural structures. 

-the existence, in space, of carbon chains and molecules 
as large as acetone is easier to understand in terms 
of breakup of large, inhomogeneous, kerogen-like, grains 
(formed in "dense" circumstellar shells and destroyed in 
strong IS shocks), rather than in terms of synthesis in the 
tenuous IS medium. 



3.7 Excitation 

The relatively high ratio of short to long wavelength 
UIB intensities initially led to the conclusion that the 
emitter temperature should be very high. Under the as- 
sumption of stochastic heating by absorption of a photon, 
and for the relatively small PAHs initially considered, cal- 
culations showed that such a temperature already requires 
the absorption of UV to far UV photons [53] . The steady 
increase of the absolute intensities of the UIBs with am- 
bient UV radiation field intensity (G value) over 4 to 5 
orders of magnitude has been considered as the strongest 
argument in favour of stochastic heating by UV photons 
[4]. These conclusions are in conflict with a number of 
astronomical and laboratory facts. 

l)The relative UIB intensities from reflection nebulae 
are independent of the temperature of the illuminating 
star, i.e. of the hardness of the radiation field, from 22000 
down to about 5000 K (see [5]); this is at variance with the 



sharp decrease of dust temperature with photon energy 
predicted by calculations [53]. 

2) Shortward of 10/im, the shape of UIB spectra of sev- 
eral galaxies that have been studied are strikingly similar 
for nuclear regions, active regions in the disks, as well as 
quiet regions [49]. A similar stability is observed among 
spectra of regions of hugely different G values, which are 
accompanied by large variations of star temperature and, 
hence, of the spectrum of the exciting radiation. Again, 
one should therefore expect notable non-linear variations 
of relative band intensities, which is not the case [4]. the 
authors of [49] conclude that the global relation between 
UIB intensity and UV illumination might be only indirect. 
This is corroborated by the observation of UIB emission 
from poor-UV regions, such as galaxy M31 [61]. 

3) The "small PAH-stochastic heating" paradigm is 
also confronted with the stability of the UIB wavelengths 
and widths against environmental variations, especially 
in view of attendent dehydrogenation, ionization or down- 
right destruction of the molecule [4, 5, 24, 39, 48, 50] (see 
Sec. 3.2). 

4) It appears from the previous subsections that all 
available spectral evidence weighs in favour of bulk, solid- 
state models of IS dust, as opposed to free- flying individ- 
ual molecules. Recent discussions involving PAHs also in- 
voke large "molecules", up to 1000 atoms or more [5, 24, 
35, 43]. Now, according to the former calculations, this 
leads to prohibitively low dust temperatures and unde- 
tectable UIB emission. One way out of this paradox is to 
hypothesize very special, but as yet not measured, elec- 
tronic properties: high photon absorption cross-sections, 
absorption cut-off sliding to the red as the particle size 
increases, etc. This becomes still more accrobatic if one 
is to account for UIBs detected in poor-UV environments 
(see [22, 54]). 

The coal model is confronted with a symmetrical puz- 
zle as long as it holds to the simplest excitation process: 
thermal equilibrium. This is only effective in strong radi- 
ation fields, in the immediate vicinity of hot stars, as is 
the case with post-AGB stars and PPNe (Pre-Planetary 
Nebulae); see [9]); but it is out of the question in RNe 
(Reflection Nebulae) and questionable in other instances. 
All models are therefore confronted with the task of con- 
ceiving a third excitation mechanism. 

A significant hint towards this goal is provided by 
the spatial distribution of the intensity of the 3.3 feature 
across the Orion Bar observed edge-on [55] : it forms a nar- 
row peak between the ionization front (delineated by the 
Br a and Pa lines) and the H 2 excitation peak (delineated 
by its ro- vibration transitions in the ground state). While 
the sharp drop towards the ionization front may perhaps 
be attributed to partial or total destruction of the carrier 
in the (hot) HII region, the drop on the opposite side is 
more problematic because it cannot be due to a decrease 
in the column density. Sellgren et al. [55] finally assigned 
it to the extinction of the stellar radiation by molecular 
hydrogen, assuming a density of 10 4 cm~ 3 for the latter. 
However, if this was the case (with a small PAH, as re- 
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quired by single-photon heating), one would also expect 
a change in the spectrum of the exciting radiation and, 
consequently, a change in the 3.3 emission feature across 
the Bar. What is observed, in fact, is a remarkable stabil- 
ity of wavelength, width, [3.4]/[3.3] and [continuum]/ [3.3] 
ratios. 

Now, the steep, large and opposite intensity variations 
of the ionic and molecular hydrogen spectral features ob- 
served to occur across the Bar suggest that, in between 
the ionic and molecular regions, their must be a narrow 
transition layer of atomic hydrogen, roughly coincident 
with the 3.3 feature peak. A similar situation has been 
described in NGC 7027 [56] and seems to have been ob- 
served at the edge of high-latitude cirruses [41, 65]. This 
led us to surmise that the excitation of the UIBs is linked 
with the presence of these atoms, rather than directly with 
the ionizing radiation [57]. Such a scenario (developed in 
the next section) naturally and straightforwardly explains 
the small transversal extent of the 3.3 layer, without hav- 
ing recourse to uncertain or qualitative assumptions as to 
the destruction of dust carriers towards the exciting star 
[41], or steep extinction of the UV towards the molecular 
cloud. 

4. Chemiluminescence 

Hydrogen radicals (free atoms) are known to be highly re- 
active, especially with hydrogen and carbon. When such 
a radical hits a hydrogeneted surface, it is likely to form, 
with an H atom of the surface, a H2 molecule, which will 
promptly escape, leaving a dangling bond. This catalytic 
reaction between two H atoms is the main route envisioned 
by astrochemists to explain the high rate of formation of 
molecular hydrogen in space: here, the solid surface pro- 
vides the third body which is necessary to ensure conser- 
vation of momentum. Since carbon is known for its high 
catalytic efficiency, our IS carbon dust is a good candi- 
date for the job, provided it can present a large fraction 
of dangling bonds per C atom. It must therefore be dis- 
ordered and porous, as are coals and kerogens. Note that 
the nascent H2 molecule carries away a large part of the 
energy made available by this exothermic reaction; it must 
therefore be ro-vibrationally excited and this may be the 
source of part of the H2* lines observed in PDR's and 
at the edge of molecular clouds, between ionized H and 
molecular gas (CO). 

Now, when another H atom approaches the dangling 
bond left on the dust, it is very likely to be captured and 
form a new strong C-H bond. The chemical (potential) 
energy of this bond, ^4.5 eV, is insufficient to dissoci- 
ate other bonds or excite higher electronic states; it is 
then available as ground-state vibrational energy for the 
atoms in the grain. Initially, this energy is essentially in the 
form of C-H stretch. However, it quickly spreads around, 
exciting the other characteristic (normal) modes of the 
given particle. The bond anharmonicity (non-linearity of 
restoring forces) opens channels for the available energy 
to be shared among the modes, which otherwise would be 



isolated. This process is called IVR (Internal Vibrational 
Redistribution), and is usually very short-lived, of order 1 
ps, after which the energy deposited in an isolated grain 
can only escape through IR radiation, with a time con- 
stant of order 0.1 s. 

Such infrared chemiluminescence, following molecular 
reaction in the gas phase was observed long ago [58] and 
often thereafter. Similar processes occur in an isolated 
gaseous molecule upon absorption of a vis/UV photon, as 
in the experiment of Williams and Leone [21], referred to 
above. The problem of theoretically predicting the chemi- 
luminescence spectrum is extremely difficult to solve be- 
cause of lack of detailed knowledge of the various couplings 
(force fields) between atoms in a given molecule: electro- 
static dipolc and higher order moments, van der Waals and 
other forces, hydrogen bonding, etc. This is compounded 
by the length of dynamical calculations in the case of large 
systems, such as those that are now envisioned as dust 
models. Fortunately, commercial software is now available 
which is dedicated to numerical simulation of molecular 
structure, dynamics and reactions. I have used such gear 
to collect quantitative, semi-quantitative and qualitative 
data on the structure of coal/kerogen model particles as 
large as 500 atoms, the interaction of H atoms with such 
particles and the resulting vibrational spectra [59, 60]. 

The main result of these computations, of interest for 
present purposes, is that, to a large extent, the predicted 
emission spectrum is very similar to the absorption spec- 
trum, independent of particle size and colour tempera- 
ture of the illuminating star. No immediate thermaliza- 
tion occurs. Moreover, since the only relaxation route is 
through IR radiation, the luminescence efficiency in the 
UIBs should be very high (of order 0.1 if it is defined as 
the ratio of IR energy emitted over deposited energy). 

Chemiluminescence thus solves the puzzles encoun- 
tered above with the two other excitation mechanisms, 
provided enough atomic H is available. If there is a lu- 
minous star near by a molecular cloud, it will give rise, 
at the edge of the cloud, to a PDR (PhotoDissociation R 
egion) such as the Orion Bar (referred to in Sec. 3.7), with 
contiguous ionized, atomic and molecular regions, succes- 
sively in the radial direction from the star. The quantity 
of chemically excited dust in the atomic layer, and hence 
the UIB intensities, will obviously depend on G, but the 
IR emission spectrum will not! This is all compatible with 
the roughly linear increase of the individual UIB intensi- 
ties versus G, as well as with the stability of their spectrum 
over 5 orders of magnitude of G [24]. Moreover, because 
the rates of dissociation of H2 molecules and ionization 
of H atoms depend on the flux of UV photons, the UIB 
intensities will depend, in absolute but not relative value, 
on the stellar colour temperature, which is not the same 
for all G's. This, together with the random variations of 
molecular density at the cloud edge, accounts for the ob- 
served excursions from strict linearity of intensity with G 
[24]. 

Of course, this does not exclude the coexistence of ther- 
mal equilibrium heating of grains if the ambient radiation 
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Fig. 7. The optimized structure of the building block (101 
atoms) of the reference macromolecule (493 atoms) mod- 
elled in [60]. 



field is strong enough. In the opposite case of no neigh- 
bouring luminous star, atomic H is still present, but in 
smaller quantities, if only for the diffuse galactic IS radia- 
tion field; this accounts for the UIB luminosity of nebulae 
at high Galactic latitudes, such as the Chameleon cloud 
or the Ursa Major cirrus [41]. 

A further observational support for 
this excitation model is the spatial co- 
incidence of the UIB emission with the 
"filament" of excited H2 emission at the edges of 
NGC 7027 [56] and of the pOph main cloud [65]: indeed, 
as indicated at the beginning of this Section, chcmilu- 
minescence induced by atomic hydrogen is genetically 
associated with the production of hydrogen molecules 
which carry away part of the chemical energy of the 
incoming H atom and are, therefore, ro-vibrationally 
excited. 

In application of these ideas, a chemical kinetic model 
was built [11] for the galaxy M31, whose Diffise ISM emits 
the UIBs although it is nearly devoid of far UV radiation 
[61]. A good numerical fit to the observed IR intensities 
was obtained assuming an atomic H density between 50 
and 200 cm~ 3 (in agreement with the range determined 
for the edge of a typical molecular cloud in our Galaxy) 
and a chemiluminescence efficiency of 0.3, in agreement 
with the numerical chemical simulation referred to above. 

The essential difference between this excitation mech- 
anism and stochastic heating does not lie in the energy 
deposition process (atomic H capture or photon absorp- 
tion), but in the unwarranted assumption at the core of 
the stochastic heating theory, that thermalization of the 
absorbed energy occurs before radiative relaxation, only be- 
cause the latter may require up to 1 s. Common wisdom 
mixes up IVR with thermalization and the last decades 
of the past century were necessary to discriminate be- 
tween various types of chaos and between various degrees 
of randomization (see bibliography in [60]). In momentum- 
configuration phase space, thermal (statistical) equilib- 
rium is the extreme case where all points of the energeti- 
cally available phase space are equally probable and ulti- 
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Fig. 8. The spectrum of the fragment of dust grain mod- 
elled in Fig. 7, 12.2 ps after capture of an H atom; 
smoothed to a resolution of 13 cm -1 . 



mately visited by the representative point of the molecular 
system. This stage is only ideal and, in fact, very difficult 
to reach on a microscopic scale and in the finite time of 
experiments, but it is easier to approach when dealing 
with macroscopic assemblies of identical particles inter- 
acting with each other and with a "thermal bath" (itself 
an ideal concept!). By contrast, experiments resolved in 
time, space and wavelength have shown that the vibra- 
tions of atoms in an isolated system may retain a high 
degree of coherence for a long time. That means that only 
parts of phase space are visited, that randomization is not 
complete and that a temperature cannot be defined. I have 
shown numerically that this is the case for a disordered 
system, even as large as 500 atoms [60]. For this purpose, 
I used a structure inspired by chemical models of coal 
[62]. A small fragment of this is shown in fig. 7, and fig. 8 
displays the frequency spectrum of its elecric dipole fluc- 
tuations after deposition of ^4.5 eV in one C-H bond. If 
the same energy were thermalized, the peak temperature 
reached by the particle would be ~ 300K and the ratio 
[3.3]/[11.3] would fall below 1% ! 

Experiments such as that of Williams and Leone [21] 
and Wagner ct al. [15] seem to show that randomization 
in PAHs is also slow after absorption of a UV photon, 
so that IR p/iotoluminescence is observed. This process 
could therefore coexist with chemiluminescence. However, 
it probably suffers from the undesirable loss of part of the 
deposited energy to vis/UV fluorescence and electronic 
continuum emission; neither of these affects chemilumines- 
cence since the chemical excitation considered here leaves 
the system in the electronic ground state. 



5. Conclusion 

One is overwhelmed by the immense diversity of both UIB 
and coal/kerogen IR spectra. Some relief is provided by 
1) the sparsness of the prominent bands, 2) their relative 



R. Papoular: On the UIBs and their excitation mechanisms 



13 



stability in wavelength, 3) the encouraging fits attempted 
between observations and models over most of the UIBs 
and extended to the UV/vis and far IR ranges, 4) the 
vast leeway for model tailoring, offered by the disordered 
structure of coals and kerogens, as opposed to the very 
restrictive aromaticity of PAHs, 5) the credibility of the 
model material, because it is so easily found in nature and 
essentially made of only the four most abundant elements 
in space (after He), and 6) the fact that this material has 
already been so thoroughly studied in the laboratory that 
even its diversity is mostly understood. 

The issue of the excitation mechanism is intimately 
linked to the issue of modelling the observed spectra. By 
invoking chemiluminescence on large particles, one avoids 
high dust temperatures and the attendent puzzles of vary- 
ing band wavelengths and widths which are not observed 
in the sky. 

Among the issues that seem now to deserve early in- 
vestigation are 

- the assignment of the "lone" 11.3-/im feature (as op- 
posed to the mono-, duo-, etc. out-of-plane C-H bends); 

- the assignment of the 12.7 /jm feature; 

-the decreasing contrast of UIBs with increasing dis- 
tance from central star, in certain objects (see [42]); 

-the experimental demonstration of IR photo- and/or 
chemiluminescence on large isolated fragments of coal or 
kerogen and measurement of their efficiencies; 

-the need for a library of complete celestial IR spectra 
(<~2 to ~16 /Ltm), each associated with a known environ- 
ment, differing from the others by the intensity of the 
radiation field or the nature of the object : RN, PPN, PN, 
PDR, GC, etc.; this would help understand, by analogy 
with coals and kerogens, the evolution of IS dust from in- 
ception in the envelopes of AGB stars to its consumption 
in protostars; 

-the measurement of the visible absorption spectrum 
of kerogen, to be compared with the DIBs. 

The coal/kerogen model can set a frame for the discus- 
sion of these and other issues. However, keeping in mind 
that these terrestrial materials were forged under special 
physical and chemical conditions, the model will no doubt 
have to be amended according to observations, as the de- 
bate goes on. 
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